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1. Non equilibrium thermodynamics in its conventional form! is characterized 
AN : al axl 
by the linearity of the phenomenological relations between fluxes ab and affi- 


nities &, 


da" : : 
2 = €, (vu, v=1,..., 2, Summation convention). (1) 


dt 


Further the connection between thermodynamic coordinates «” and affinities is 
linear 


Si — Sip ® (8p — —— 72> (&% = 0); S(x)=entropy function). (2) 
Oa 


The time reversal invariance of the fundamental mechanical equations of motion 
is reflected by Onsager’s reciprocal relations for the phenomenological coeffi- 
cients LL” 
Ti, Ty TAY = 1 (Th, Te = 8h), (3) 
where 7! is the representative of the time reversal operation, that is with a 
suitable choice of coordinate system 
+1 (x even) 


TY =0 (ur), T= : (4) 
—] (x odd) 


Sometimes the linear approximation (1), (2) is invalid and in the pheno- 
menological equations 


a Am) (5) 


the functions a(«) are not linear. It is natural to ask whether in the case (5) 
of non linear relaxation the time reversal invariance of the mechanical equa- 
tions of motion will give rise to generalized reciprocal relations which reduce to 


1 B.g. [de Groot, 1951]. 
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Onsager’s relations in the linear case. In this note we shall show that as a 
consequence of the time reversal invariance the equation (5) can be written 
0 log f(x) , aL’ (a) 


» LLM (a) f (a2)] = LE (a) a (6) 


dal 4) 
te 
(ollod 


where f(x) is the density of the thermodynamic equilibrium distribution of the 
g-observables and L'”(«) satisfies 


THT. EF (Pa) = L™ (a), (7) 
where 7! has the same meaning as above. The equilibrium distribution is in- 
o(T «) 
variant under time reversal and as the Jacobian 3 (a) =, we obtain 
a 
f (Da) =f (x). (8) 


If in particular L(x) is independent of « and if f(«) has the Gaussian form 
f (a) ~ exp { — $ Spy a! a” }? (9) 


it is seen that (6) is identical with (2), and (7) with (3). 


2. The proof of (6) belongs to statistical mechanics. We have associated with 
the g-observables an equilibrium distribution—not necessarily a probability dis- 
tribution— 


fla)da (10) 


in the n-dimensional real space of the thermodynamic coordinates. The tempo- 
ral change of the «-observables is partly characterized by the transition prob- 
ability distribution 


IPS, CR, 10). (11) 


The expression (11) represents the probability of finding at the time ¢ a value 
of the observable in the n-dimensional interval J, if at the previous time s the 
value of the observable is known to be the n-dimensional point «. 

The transition probability distribution (11) is assumed to have the following 
properties (on a macroscopic time scale). 

a. For fixed real values s,t(s<t) and fixed « in the n-dimensional real space 
k", P(s,a;t,.) determines a probability distribution in R”. 

b. For fixed s, t, and J the derivatives 


a) 
ea eth 1) Rael Aa o — P(s,a;t, 1) (12) 


exist as continuous functions of «. 
c. Continuity condition. For arbitrary 6 >0 


! Boltzmann’s constant = 1. 


362 


ARKIV FOR FYSIK. Bd 17 nr 20 


1 ; 
hi IP (S,.6& )V= Oz ‘ 
\ ee { (s,ost,da)=0. (13) 
|a’-a|>3d 
It follows from this assumption that 
1 (2 eT) 
lin (Sen tel y= ; (14) 
e940 0 («¢Z) 


d. Kolmogorov’s conditions.” For arbitrary 6>0 the following limits exist 


. 1 i , , 
lim ; | Pe otd alla x)= (ot): (15) 
(t-s)}ot —S “) 
|a’—a|<d 
- 1 i rf tu 1 ‘yp y iw 
lim Piero. da) (o" ae) (a a) 10 \(o,,0)) (16) 
C= tOt—Ss A 
|@’=a)|<0 


It follows from the continuity condition c. that the limits a(«,t) and b!” («,t) 
are in fact independent of 6. Further it follows that for arbitrary real numbers 
k,,...,kn the expression 
BP (a,b) Mi ey (17) 
is non negative. 
Let g(.) be an arbitrary complex function on R”, continuous and bounded 


o a2 

; : aes Cp (a) a py (a) 

and having continuous derivatives ao and - - 
Ca On Ca 


,- The conditioning transfor- 
mation C(s,t) is defined by 


C(s,t)@ (a)=[P(s, «5 t,da’')p(a’), (18) 


R 


and represents the conditional expectation of the function @(.) of the observable 
at the time ¢ given the value « of the observable at the time s. Evidently the 
transform C(s,t)p of m is according to the assumption b above bounded and 
twice continuously differentiable. 

It follows by standard methods? from the continuity condition and from 
Kolmogoroy’s conditions that the limit 


é ee ACN ae a 
[Seen] = tim, ae p (a) 


~ 


== han el P(s,at.da)\[p(a)—m(a)] (19) 
(t-s)yot—Ss 
RY 


=n 
2 pe |2 
+ faP= 2 lat. 
= . . 
2 E.g. [Gnedenko, 1958]. Kolmogorov’s conditions are usually associated with Markov pro- 


cesses of the s.c. continuous type. Here, however we are not assuming that the transition prob- 
ability distribution (11) can correspond to a Markov process. 
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exists and is equal to 


Oy (x) ' O° ep (a) 9 
= Ait Lb (a. ty —— 2) 
M (t) @ (a) =a" (a, t) ah 4b (a, ) a age (20) 
or in operator form 
(Sy te) — Ul 
ene OM aes (21) 
(t-s)y0 =a 


In the following we shall call M(t) the transport operator corresponding to the 
a-observables. 

e. We may restrict ourselves to the temporally homogeneous case so that 
the transition probability distribution P(s,x;t,/) is a function only of the time 
difference (t—s) 


Psat) =P(ast—s-f), (22) 
or equivalently 
C(s,t)=C(t—s). (23) 
In this case the transport operator M(.) is independent of the time 


M(t)=M. (24) 


The following transformation property of the transport operator M is some- 
times useful. Let 


w=ul(a),...% =u" (a) (25) 


be m functions of the n-dimensional variable x, having continuous partial de- 
rivatives of the second order. Applying the transport operator M to the function 


p(x) =p (ut (a), ..., &” (a) (26) 
we obtain directly by formal calculations 


du! du" ey 
Oa“da’ awau* (27) 


OW 
Mo= (Mw) +0" 
ACI) 


3. We have assumed that the distribution (10) corresponds to equilibrium. 
This imphes that the transition probability distribution (11) admits the distri- 
bution (10) as a stationary distribution, that is 


[i (a)daP (at, 1)=[fla)da. (28) 


R” I 
Hence for an “arbitrary” function @ («) 


ja p (x) fla)da=-0. (29) 


R 
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By formal partial integration we find that (29) is equivalent to 


7) 1 
a (at) f (ce) — 5 


ra) 
role Oa 


ec) (| = 0; (30) 
It is seen at once from (30) that the expression 


J" (ce) = a (ct) f (ce) — 5 = B” (x) f (x) (31) 


may be heuristically considered as the component of a flow vector in the n- 
dimensional real space. The relation (30) then expresses the fact that the flow 
is mcompressible (divergence free). The incompressibility relation (30) is obviously 
satisfied automatically if we make the generally much stronger assumption of 
detailed equilibrium 


J” (x) =0. (32) 


It follows from (31) that in the case of detailed equilibrium the coefficient a’ (x) 
is determined by f(x) and b””’ («) 


ma 


al (2) =f (2) Soe (AY (2) (2). ce) 


This relation is indeed of the form (5), (6). We shall see that this relation can 
really be derived from statistical mechanics. It may be noticed however that 
in the one-dimensional case (7=1) detailed equilibrium holds already if, as is 
frequently the case in the applications 


[fla|a(a)|da< co, lim f(a)b(a)=0. (34) 
|ale~ 


—o 


For, integrating (30) we find 


a(x) f («) -~-~— f (a) b (a) = const. (35) 


Integrating both sides of the equation over an unbounded interval we obtain 
on the left a finite value whereas on the right we obtain a finite value if and 
only if the constant is chosen to be zero. Hence (32) must be satisfied. 


4. For the further study of the transport operator M it is convenient to in- 
troduce the Hilbert space L* of complex functions of the n-dimensional real 
variable «, X(«), Y(«),... square integrable with respect to the measure f(«)d«. 
The scalar product in L? is defined by 


(X| Y>=[X (a) V (a) f (a) da (36) 


R 


365 


U. UHLHORN, Onsager’s reciprocal relations for non linear systems 


and the norm of an element X of L? by 


|| X ]=< x 


Cus (37) 


From the invariance of f(x)d«, (28) it follows at once by application of Schwarz’ 
inequality that the conditioning operator C(t) is bounded by one 


\| C(t) X || <|| XI]. (38) 


As a consequence of the invariance (8) the time reversal transformation 7’ in 
I? defined by 


TX (a) =X (T «) (39) 
is a unitary involution 
(Bes Maeaed bee (40) 
the * denoting Hermitian adjoint. 


The time reversal invariance of the mechanical equations of motion gives 
rise to the symmetry 


like )\daP(a;t,L:)={ f(a)daP (a;t,T I.) (41) 


TI, 


for arbitrary n-dimensional intervals J,;, /;; ZZ denoting the “‘time reversed inter- 
val’. In terms of the time reversal transformation (39) and the conditioning 
operator C(t) the relation (41) can be expressed as 


TC (t)T =C(-t)=C (t)*. (42) 

Hence we find that the transport operator M_ satisfies 
TMT =M*. (43) 
The relation (43) is symbolic as the operator M, being unbounded, is not 


defined in the whole of L*. Howewer its domain of definition is dense in L?. 
We now procede to determine M*. By the definition and using (27) 


=<(MX|Y)>={ (MX) Yfde= 


R” 


=|(MXY)fd«e-[X(M VY) fda— [oe 
R” Ry R” 
Hence 


0 0 
M*ty= M 14 ae bi” eng ype 
/ <a| tp 


¥ eye Ox xi BY. 
-|- a! +f J otwep |e } . bi ra) } 
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In order to obtain 7 MT we consider 


TMT X (x) =MTX (Ta) =a" (Tx) = L pi pig) 2 EX 
(x) AT X(T «) =a" (T «) Aol (Ta) + 5b (Do), mag %) 

Posen 1 ox \ l uy" t v aX 

Cd, (T a) Die 5 yn (&) + 5 (L a) Dy Py aoe (46) 
Equating (45) and (46) we find 
Ve 1 4 = C m 

a (La) The = — al (x) + f(a) * 5 (O(a) F(a] (47) 
BE (Ligh) Tele = Oo) (48) 


We recall that equation (30) has the form of an incompressibility relation 


dd (a) 
—- =(0,; (49) 


18 (cl (ax) f ( 
i ee (50) 
where 
C" (a= —c" (a)-. (51) 


Hence we find the following representation for a" («) 


f(a) a (a) =~ LL!" (2) f(a) (52) 
where 
TM (x) = 4 [0 (x) +e” (a). ee) 


From (52) we obtain easily the following representations of M and M™*, (45): 


; alee : 

ux=p'%, fz 2 | (54) 
0 0X. 

M*X = cna [ Le a . (55) 


We propose that it is possible to choose ¢"(«) so that L'’(«) satisfies the 
symmetry relation (7). As is seen from (48) the symmetry relation is valid for 
the symmetric part of L/”(x), ie. $6" («). It is therefore sufficient to show that 
c’” (x) can be chosen so that it satisfies (7). Assume that c,"” (%) satisfies (50) 


and define c’”(«) by 


1 This is wellknown i the 3-dimensional case: If V- ee 0 we can write ne V x Ae 
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co” (x) = 4 [er (a) +r” (Lax) Th Ty). (56) 


It is evident that c!’(x) defined by (56) satisfies the symmetry relation (7) so 
that we have only to prove that it satisfies (50). We find, using (47), after an 
elementary calculation 


es 


2 Uh (a) oF (Po) Pi Te] = 5 (a) (a). (57) 


20 


and this completes the proof. 


5. The derivation of the reciprocal relations (7) for the non linear system (6) 
has been based on the time reversal invariance of the mechanical equations of 
motion and on the assumptions about the behaviour of the transition proba- 
bility distribution (11) for macroscopically small time differences expressed by 
Kolmogorov’s condition. It should be noticed that we have made no assumptions 
about the behaviour of the transition probability distribution for large time 
differences. In particular we have not assumed that the transition probability 
distribution corresponds to a Markov process. It is evident that the assumptions 
underlying Kolmogorov’s second condition are compatible with statistical mechan- 
ics only on a sufficiently coarse grained time scale. This is evidently the same 
assumption as used in the linear theory (see in particular [Smoluchowski, 1916]). 
The conditions of Kolmogorov determine the behaviour of the transition prob- 
ability distribution for macroscopically small time differences only, in fact they 
determine only the derivative with respect to time at the origin of time, whereas 
all further properties of the transition probability distribution remain unspeci- 
fied and are irrelevant to our discussion. 


ACKNOWLEDGEMENT 


IT wish to thank Professor Oskar Klein, head of this institute, for his encouragement and 
unfailing interest in the work. 


Institute for theoretical physics, University of Stockholm, Stockholm. February 1960. 


REFERENCES 


GNEDENKO, B. W., 1958, Lehrbuch der Wahrscheinlichkeitsrechnung, Akademie-Verlag. 

bE Groot, S. R., 1951, Thermodynamics of irreversible processes, North Holland Publishing 
Company. 

SMoLucHOwSsKI, M., RirreR von Smoxav, 1916, Drei Vortraige tiber Diffusion, Brownsche Mole- 
kular bewegung und Koagulation von Kolloidteilchen, Phys. Zeitschr. 17 557, 585. 


Tryckt den 31 mars 1960 


Uppsala 1960. Almqvist & Wiksells Boktryckeri AB 


